ITP-Budapest Report No. 635 



^ : Form factors in finite volume II: 

^ ; disconnected terms and finite temperature correlators 

^ ■ B. Pozsgay^* and G. Takacs^^ 



in 



o 



X 



^Institute for Theoretical Physics 
Eotvos University, Budapest 



"^HAS Research Group for Theoretical Physics 
'.st, Pdzmdny Pe 

25th June 2007 



^! H-1117 Budapest, Pdzmdny Peter setdny 1/A 



Abstract 



> 

in 
o 

^ , 

\ Continuing the investigation started in a previous work, we consider form factors of 

\^ • integrable quantum field theories in finite volume, extending our investigation to matrix 

\ elements with disconnected pieces. Numerical verification of our results is provided by 

■ truncated conformal space approach. Such matrix elements are important in computing 

finite temperature correlation functions, and we give a new method for generating a low 
temperature expansion, which we test for the one-point function up to third order. 



■ 1 Introduction 



The matrix elements of local operators, the so-called form factors are central objects in quan- 
tum field theory. In two-dimensional integrable quantum field theory, the S matrix can be 
obtained exactly in the framework of factorized scattering (see [1, 2] for reviews). Using the 
scattering amplitudes as input, it is possible to obtain a set of axioms [3] which provides the 
basis for the form factor bootstrap (see [4] for a review). 

Although in the bootstrap approach the connection with the Lagrangian formulation of 
quantum field theory is rather indirect, it is thought that the general solution of the form 
factor axioms determines the complete local operator algebra of the theory [5], which was 
confirmed in many cases by explicit comparison of the space of solutions to the spectrum of 
local operators [6, 7, 8, 9]. Another important piece of information comes from correlation 
functions: using form factors, a spectral representation for the correlation functions can be 
built which provides a large distance expansion [10, 11], while the Lagrangian or perturbed 
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conformal field theory formulation allows one to obtain a short-distance expansion, which 
can then be compared provided there is an overlap between their regimes of vaUdity [11]. 
Other evidence for the correspondence between the field theory and the solutions of the form 
factor bootstrap results from evaluating sum rules like Zamolodchikov's c-theorem [12, 13] or 
the A-theorem [14], both of which can be used to express conformal data as spectral sums 
in terms of form-factors. Direct comparisons with multi-particle matrix elements are not so 
readily available, except for perturbative or 1/N calculations in some simple cases [3]. One 
of our aims is to provide non-perturbative evaluation of form factors from the Hamiltonian 
formulation, which then allows for a direct comparison with solutions of the form factor axioms. 

Based on what we learned from our previous investigation of decay rates in finite volume 
[15], in our previous paper [16] we determined form factors using a formulation of tlie field the- 
ory in finite volume. We used the truncated conformal space approach (TCSA) developed by 
Yurov and Al.B. Zamolodchikov [17] as a basis for numerical comparison to non-perturbative 
Hamiltonian formulation of quantum field theory, and also its fermionic version in the case 
of the Ising model [18]. We were able to give an extensive and direct numerical comparison 
between bootstrap results for form factors and matrix elements evaluated non-perturbatively. 
One of the advantages is that we can compare matrix elements directly, without using any 
proxy (such as a two-point function or a sum rule); the other is the very high precision of the 
comparison and also that it is possible to test form factors of many particles which have never 
been tested using spectral sums. Our approach, in contrast, makes it possible to test entire 
one-dimensional sections of the form factor functions using the volume as a parameter, and 
the number of available sections only depends on our ability to identify multi-particle states 
in finite volume. Part of the motivation of this work is to complete the non-perturbative eval- 
uation of form factors by extending our results to matrix elements with disconnected pieces. 

Another motivation is provided by the fact that such matrix elements are relevant for 
the calculation of finite temperature correlators. Finite temperature correlation functions 
have attracted quite a lot of interest recently [19, 20, 21, 22, 23, 24, 25, 26, 27]. Leclair 
and Mussardo proposed an expansion for the one-point and two-point functions in terms of 
form factors dressed by appropriate occupation number factors containing the pseudo-energy 
function from the thermodynamical Bethe Ansatz [20]. It was shown by Saleur [21] that 
their proposal for the two-point function is incorrect; on the other hand, he gave a proof of 
the Leclair-Mussardo formula for one-point functions provided the operator considered is the 
density of some local conserved charge. His proof is based on a conjecture concerning the 
expression of diagonal finite volume matrix elements in terms of connected form factors. In 
view of the evidence it is now generally accepted that the conjecture made by Leclair and 
Mussardo for tlie one-point functions is correct; in contrast, the case of two-point functions 
(and also higher ones) is not yet fully understood (see the introductory part of section 7 for 
more details). Here we investigate how finite temperature one-point functions can be expanded 
systematically using finite volume L as a regulator and make a proposal which is expected to 
be valid for multi-point correlators as well. 

Our exposition is structured as follows. In section 2, after recalling the form factor boot- 
strap axioms, we present a brief review of the approach developed in our earlier paper [16] 
(to which we refer the interested reader for more details), and then we state our main re- 
sults which is the description of all matrix elements containing disconnected contributions. 
In Section 3 we briefly recall the two models used for numerical comparison, which are the 
scaling Lee- Yang model and the Ising model in a magnetic field. We omit the description of 
the method for obtaining matrix elements from truncated conformal space, and instead we 
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refer the interested reader to [16] where all the necessary details can be found. 

As we showed in [16], there are essentially two types of matrix elements with disconnected 
contributions. Section 4 is devoted to the first type, which is the case of diagonal matrix 
elements; we present a general formula for them in terms of the symmetric evaluation of the 
diagonal form factor and test it against truncated conformal space. In section 5 we analyze 
diagonal matrix elements in terms of connected form factor amplitudes, and we show that 
our results are fully consistent with the above-mentioned conjecture made by Saleur in [21]. 
In section 6 we discuss the second type of matrix elements with disconnected contributions, 
namely those with particles of exactly zero momentum in the finite volume states. Adding the 
results presented in section 4 and section 6 to those obtained in [16], we achieve a complete 
description of all multi-particle matrix elements of a general local operator to all orders in 
1/L. Section 7 is devoted to finite temperature correlation functions: we propose a systematic 
method for deriving a low-temperature expansion, which is applied to one-point functions and 
tested by comparing the results to the Leclair-Mussardo expansion [20]. We also briefly discuss 
the extension of our method to the evaluation of two-point functions. Section 8 is reserved for 
the conclusions. 



2 Form factors in finite volume: a brief review 
2.1 Form factor bootstrap 

Here we give a very brief summary of the axioms of the form factor bootstrap, because we 
need them in the sequel; for more details we refer to Smirnov's review [4]. Let us suppose for 
simplicity that the theory has particles Ai, i = 1,. . . ,N with masses which are strictly 
non-degenerate i.e. rrii / rrij for any i j (and therefore the particles are also self-conjugate). 
Because of integrabihty, multi-particle scattering amphtudes factorize into the product of 
pairwise two-particle scatterings, which are purely elastic (in other words: diagonal). This 
means that any two-particle scattering amplitude is a pure phase, which we denote by Sij (6) 
where 9 is the relative rapidity of the incoming particles and Aj. Incoming and outgoing 
asymptotic states can be distinguished by the ordering of the rapidities: 

m a\ _ j\^l^---^^n)il..i„ : 9i > 62 > ■ ■ ■ > 0n 

' ll^l,---,^n)r.i„ : ^l<^2<---<en 

and states which only differ in the order of rapidities are related by 
1^1, ...,9k, 9k+l, . . . , dn)ii...ikik+i-in = Si^,i^j^^{9k — 0fe+l)|0l, . . . , 9k+l,9k, ■ ■ ■ , 9n)n...ik+iik--in 

The normalization of these states is specified by giving the following inner product among 
one-particle state: 

j{9'\9)i = 6ij27Td{e' -9) 
For a local operator 0{t,x) the form factors are defined as 

-^mn(^m> ■ ■ ■ >^ll^l>- • • > ^n)ji...j,„;ii...i„ = ji-jmi^l, ■ ■ ^mP(0' 0)1^1' • • ■ ,0n)ii...in (2-1) 
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With the help of the crossing relations 

-^mn(^l' • • • '^ml^l' • • • } ^n)ji...jm;ii---in ~ 

-^m-ln+l(^l) • • • ) ^m-ll^m + ^'^i ^1) • • • ) Sri)ji...jm-r,3mh---in 
n k—1 

+ E "^^^jmik^ie'^ - Ok) n Sni, {01 - Ok) 

k=l 1=1 
^-^m-ln-l(^l! • • • '^m-ll^l) • • • )^A:-1,^A:+1 • • • i0n)h...jm-V,3mn-ik-iik+i-in (2-2) 

all form factors can be expressed in terms of the elementary form factors 

F^{eu 9n)n...in = (0P(0, 0)1^1, ... , en)n...ir. (2.3) 

which satisfy the following axioms: 

I. Exchange: 

{Oil ■ ■ ■ -i^ki Ok+l-i ■ ■ ■ ) ^n)n...ifcifc+i...in ~ 
Sikik+ii^k - Ok+l)F°{ei, . . .,0k+l,0k, . . ■ ,On)h...ik+iik-in (2-4) 

II. Cyclic permutation: 

(01 + 2i7r, 02,..., On) = (02, ...,en, 9,) (2.5) 

III. Kinematical singularity 

-i ResFf+2(^ + ■ ■ ■,On)ijh...i„ = ( 1 - ^ij fl ^^^ki^ - Ok)] i^f (^l, • • ■,On)h...i„ 

V k=i J 

(2.6) 

IV. Dynamical singularity 

~i Res F^^^ie + iuy2, e' - iu\j2, 0i, . . . , 0„)i,-,,.,„ = rf,if+i(0, 01, ... , 0n)fcn...i„ (2-7) 

whenever k occurs as the bound state of the particles i and _7, corresponding to a bound state 
pole of the S matrix of the form 

where is the on-shell three-particle coupling and it^^- is the so-called fusion angle. The 
fusion angles satisfy 

m| = -|- -|- 2mimj cos ufj 
27r = ufj + uj^ + uik 



and we also used the notation itf^- = tt — u^j. The axioms I-IV are supplemented by the 
assumption of maximum analyticity (i.e. that the form factors are meromorphic functions 
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which only have the singularities prescribed by the axioms) and possible further conditions 
expressing properties of the particular operator whose form factors are sought. 

We remark that with the exception of free bosonic theories, all known exact S matrices 
satisfy 

SiiiO) = -1 

and therefore the elementary form factors (2.3) have an exclusion property: they vanish when- 
ever the rapidities of two particles belonging to the same species coincide. 

2.2 Finite volume matrix elements to all orders in 1/L 

Following our conventions in [16], the finite volume multi-particle states can be denoted 

. . . , In})h...in,L 

where the Ik are momentum quantum numbers and ik are particle species labels. We order 
the momentum quantum numbers in a monotonically decreasing sequence: /„ > ••• > Ii, 
which is just a matter of convention. The corresponding energy levels are determined by the 
Bethe-Yang equations 

Qk{0i,...,0n) = mii^Lsmhek + ^6i^ii{ek-9i) = 2TrIk , k = l,...,n (2.9) 

l^k 

which must be solved with respect to the particle rapidities Ok, where 

Sij{e) = -Hog Sij{e) 

are the two-particle scattering phase-shifts and the energy (with respect to the finite volume 
vacuum state) can be computed as 



rriik cosh 4 



fc=i 

The density of n-particle states can be calculated as 

Pn..,AOi,...,0n)=detj(-^ , jin) ^ dQk(e^_^ ^ k,l = l,...,n (2.10) 

We are interested in matrix elements of local operators between finite volume multi-particle 
states: 

n...jj{l'l, ■ ■ ■ , l'm}\O{0, 0)|{/i, . . . , /„})i,..,„,L 

which can be obtained numerically using truncated conformal space (for details see [16], section 
3.3). On the other hand, using our previous results (eqn. (2.16) of [16]), the finite volume 
behaviour of local matrix elements can also be given as 



Pii...i„{Ol, ■ ■ ■ Jn)Pji...jmi^l^ • • • '^m) 



+ 0{e-'''^) (2.11) 
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and 9k {d'j.) are the solutions of the Bethe-Yang equations (2.9) corresponding to the state 
with the specified quantum numbers /i, . . . , (/{,..., /'J at the given vohnne L. The above 
relation is valid provided there are no disconnected terms i.e. the left and the right states do 
not contain particles with the same species and rapidity: the sets |(ii,6'i), . . . , {imOn)^ and 

{(ji, l^i); • • • > {jm, 6*^)} are disjoint. 

We recall from [16] that eqns. (2.9,2.11) are exact to all orders of powers in 1/L; we refer 
to the corrections non-analytic in 1/L (eventually, as indicated, decaying exponentially) as 
residual finite size effects, following the terminology introduced in [15]. 

2.3 Disconnected contributions 

Let us consider a matrix element of the form 



.(K,...,0|o(o,o)|{7i,...,/4),,..,„,L 



Disconnected terms appear when there is at least one particle in the state on the left which 
occurs in the state on the right with exactly the same rapidity. The rapidities of particles as 
a function of the volume are determined by the Bethe-Yang equations (2.9) 



and 



Qk{h,---,0n)=muLsmh.ek + ^5i^ii{0k-0i) = 2'Klk , k = l, 

l+k 

Qk{e[,...J'^)=mj,Lsmhei + Y,Sj.jM-0[) = 27rli , k = 1, 



, m 



Due to the presence of the interaction terms containing the phase shift functions S, equality of 
two quantum numbers 1^ and does not mean that the two rapidities themselves are equal 
in finite volume L. It is easy to see that in the presence of nontrivial scattering there are only 
two cases when exact equality of the rapidities can occur: 

1. The two states are identical, i.e. n = m and 

{jl---jm} = {il...in} 

In section 4 we show that the corresponding diagonal matrix element can be written as 
a sum over all bipartite divisions of the set of the n particles involved (including the 
trivial ones when A is the empty set or the complete set {1, . . . ,n}) 

^,..,A{h...In]\0\{h...In])^,...^^,L ^ ^ 



p({l,... ,n})L 

HA)LP{{l,...M\A)L + 0{e-i^^) 

Ac{l,2,...n} 

where |^| denotes the cardinal number (number of elements) of the set A 

p{{ki, ...,kr})L= Pik^...ik^{Oki,- --Jkr) 
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is the r-particle Bethe-Yang Jacobi determinant (2.10) involving only the r-element 
subset 1 < ki < ■ ■ ■ < kr < n oi the n particles, and 



Fiiiei, ei)i,...i, = lim F§{9i + ZTT + e, . . . , 01 + ZTT + e, 6*1, . . . , 0/)^ 



J^{{ki,...,kr})L = F2ri^ki,---,0kr)iki-ikr 

a. _i_ ^--7^ _i_ ^ a. 0. , 

lil-.-ilH-.-il 

is the so-called symmetric evaluation of diagonal multi-particle matrix elements. 
2. Both states are parity symmetric states in the spin zero sector, i.e. 

{IlT--,In} = {— /n, ■ • • 5 — A} 
{I'll ■ ■ ■ ^I'm} = {~/m, ■ ■ ■ , "/(} 

and the particle species labels are also compatible with the symmetry, i.e. in+i-r = V 
and jjn+i-r = jr- Furthermore, both states must contain one (or possibly more, in a 
theory with more than one species) particle of quantum number 0, whose rapidity is then 
exactly for any value of the volume L due to the symmetric assignment of quantum 
numbers. In section 5 we state the following conjecture 

/2fc+i,2/+i = ({/( , ...,4,0, -4, . . . , -4}|$|{/i, ...,ii,o, -ii, . . . , -h})L 

1 



P2k+i{0[, ... ,(9^,0, ... , -9[)p2i+i{0i, ... ,01,0, -9i, ... , -9i) 

[j'kA^'i^ ■ ■ ■ , ^k\^i, ■■■A)+mL F2fc+2z(ivr + 9[,...,m + 9^ 

m-9'k,...,m-~9[A,---, k -k -Oi)) + ©(e"/^^) 

where p„ is a shorthand notation for the n-particle Bethe-Yang density (2.10) and equal- 
ity is understood up to phase conventions (cf. section 5) and 

^k,liGl, ■ ■ ■ ■ ■ ■ = 

1^ P2k+2i+2i^^ + 9[ + e, . . . ,iTr + 9k + e,iTr - 9k + e, . . . ,iTr - 9[ + e, 
in + €,0,9i,...,9i,-9i,...,-ei) 

is defined by assigning the same shift e to all rapidities entering the left (or equivalently 
the right) state and taking the limit e — > 0. For the sake of simplicity we assumed 
above that there is a single particle species with mass m, but the prescription can be 
easily extended to theories with more than one particle species; an example is shown in 
subsection 7.2. 

3 Exact form factors 
3.1 Scaling Lee- Yang model 

The Hamiltonian of scaling Lee- Yang model takes the following form in the perturbed confor- 
mal field theory framework: 

L 



Jo 
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where 




is the conformal Hamiltonian and $ is the only nontrivial primary field, which has conformal 
weights A = A = —1/5. When A > the theory above has a single particle in its spectrum 
with mass m that can be related to the coupling constant as [28] 

A = 0.09704845636- • • x m^^/^ 



and the bulk energy density is given by 



The ^-matrix reads [29] 



B = -§m' (3.1) 

SMO) = ""^'^^""1 (3.2) 
sinh ^ — z sin ^ 

and the particle occurs as a bound state of itself at ^ = 2iri/3 with the three-particle coupling 
given by 

= -2^3 

where the negative sign is due to the nonunitarity of the model. In this model we define the 
phase-shift via the relation 

SLY{e) = -e^^(^) 

so that 6{0) = 0. This means a redefinition of Bethe quantum numbers Ik in the Bethe-Yang 
equations (2.10) such they become half-integers for states composed of an even number of 
particles; it also means that in the large volume limit, particle momenta become 

• UP 27r4 
m sinh tfk = — ? — 
ij 

Form factors of the trace of the stress-energy tensor were computed by Al.B. Zamolodchikov 
in [11], and using the relation 

e = zA7r(l - A)$ 



we can rewrite them in terms of They have the form 

Fni9i,...,9n) = {^)HnQnixi,...,Xn)f[ f[ l^^l^^ 

i=i j=i+i ■> 



(3.3) 



with the notations 



cosh 6* — 1 ,. „x /. 
V Jo tsinh irt 
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and the exact vacuum expectation value of the field $ is 

($) = 1.239394325 ■■■ xi m"^/^ 

The functions Q„ are symmetric polynomials in the variables Xj. Defining the elementary 
symmetric polynomials of n variables by the relations 



J](x + Xi) = ^ x"- Vf ^ (xi , . . . , Xn) , ) = for z > n 



i=l i=0 

they can be constructed as 



,n 1 /I (2) (3) (3) 

Qn = ^i"^^n-lPn , n>3 

Pn = detM^"-^ where M\f = 4"-2j+i > i, j = 1, . . . , n - 3 
3.2 Ising model with magnetic perturbation 

The critical Ising model is the described by the conformal field theory with c = 1/2 and lias 
two nontrivial primary fields: the spin operator a with A^- = A^- = 1/16 and the energy 
density e with Ag = A^ = 1/2. The magnetic perturbation, defined using the Hamiltonian 
(where Hq denotes the Hamiltonian of the c = 1/2 conformal field theory) 



H = H^ + h [ dxa{0, x) 

^0 



is massive (and its physics does not depend on the sign of the external magnetic field h). The 
spectrum and the exact S matrix is described by the famous Eg factorized scattering theory 
[30], which contains eight particles Ai, i = 1, . . . ,8 with known mass ratios, and the mass gap 
relation is [31] 

mi = (4.40490857... )|/i|^/^5 



or 



h = Khml^^^ , Kfe = 0.06203236 .. . (3.4) 
The bulk energy density is given by 

B = -0.06172858982 • • • x (3.5) 

We also quote the scattering phase shift of two Ai particles for A = 0, which has the form 

fll fll f2l sinh + z sin TTX , 

All the other amplitudes Sab are determined by the S matrix bootstrap [30]; we only quote 
the Ai — A2 scattering amplitude 



Si2i9) 



5 f J 15 (J 5 (J 15 
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because it enters some matrix elements examined later. In this model we define the phase- 
shifts by the relations (for detailed explanation cf. [16]) 

5n(^) = -e^'^"(^) and Si2(0) = e^''^^^^^ 

so that again (5ii(0) = 612(0) = 0. The form factors of the operator e in the £^8 model were 
first calculated in [32] and their determination was carried further in [33]. The exact vacuum 
expectation value of the field e is given by [34] 

{e) = eh\hf^' , 6^ = 2.00314... 

or in terms of the mass scale m = mi 

(e) = 0.45475 • • • x m 

For practical evaluation of form factors we used the results computed by Delfino, Grinza and 
Mussardo, which can be downloaded from the Web in Mathematica format [35]. They use the 
following normalized operator: 

and so all data we plot in the sequel are understood with the same normalization. 



4 Diagonal matrix elements 

4.1 Form factor perturbation theory and disconnected contributions 

In the framework of conformal perturbation theory, we consider a model with the action 

Aifi, \) = AcFT - fJ- J dtdx^{t,x) - X J dtdx^{t,x) (4.1) 

such that in the absence of the coupling A, the model defined by the action A{ii, A = 0) is 
integrable. The two perturbing fields are taken as scaling fields of the ultraviolet limiting 
conformal field theory, with left/right conformal weights /i$ = /i$ < 1 and /7.<i> = /7.<i> < 1, i.e. 
they are relevant and have zero conformal spin, resulting in a Lorentz-invariant field theory. 

The integrable limit A{ii, A = 0) is supposed to define a massive spectrum, with the 
scale set by the dimensionful coupling /j,. The exact spectrum in this case consists of some 
massive particles, forming a factorized scattering theory with known S matrix amplitudes, 
and characterized by a mass scale M (which we take as the mass of the fundamental particle 
generating the bootstrap), which is related to the coupling /x via the mass gap relation 

where k is a (non-perturbative) dimensionless constant. 

Switching on a second independent coupling A in general spoils integrability, deforms the 
mass spectrum and the S matrix, and in particular allows decay of the particles which are 
stable at the integrable point. One way to approach the dynamics of the model is the form 
factor perturbation theory proposed in [36]. Let us denote the form factors of the operator ^' 
in the A = theory by 

(01, ... , ^?„),^ = (0|*(0, 0)1^1 . . . 0„)^=o„ 
Using perturbation theory to first order in A, the following quantities can be calculated [36]: 
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1. The vacuum energy density is shifted by an amount 

SS^ac = A (0| ^ |0);^=o • (4-2) 

2. The mass (squared) matrix M^^ gets a correction 

5Ml, = 2\F^{m,Q)^^,5^^,m, (4.3) 

(where the bar denotes the antiparticle) supposing that the original mass matrix was 
diagonal and of the form M^^ = m^dab ■ 

3. The scattering amplitude for the four particle process a + b — >■ c + dis modified by 

^5-^(^,A) = -^A ^'(^"'^ + ^"'y^--^-°^ , e = e.-e,. (4.4) 

marrib smh 

It is important to stress that the form factor amplitude in the above expression must be 
defined as the so-called "symmetric" evaluation 

Jim if {iTT + e,e + iTr + e,0, 0)-^^^ 

(see eqn. (4.9) below). It is also necessary to keep in mind that eqn. (4.4) gives the 
variation of the scattering phase when the center-of-mass energy (or, the Mandelstam 
variable s) is kept fixed [36]. Therefore, in terms of rapidity variables, this variation 
corresponds to the following: 



5S'J,{e,X) = ^^^ '-59 + \ 



89 d\ 

where 



A=0 



rriaSma + maSnia + {mbdnia + maSmb) cosh 9 
= m„m, sinh^ 

is the shift of the rapidity variable induced by the mass corrections given by eqn. (4.3). 

It is also possible to calculate the (partial) decay width of particles [33], but we do not need 
it here. 

We can use the above results to calculate diagonal matrix elements involving one particle. 
For simplicity we present the derivation for a theory with a single particle species. Let us start 
with the one-particle case. The variation of the energy of a stationary one-particle state with 
respect to the vacuum (i.e. the finite volume particle mass) can be expressed as the difference 
between the first order perturbative results for the one-particle and vacuum states in volume 
L: 

Am(L) = XL (({0}|*|{0})l - (0|*|0)i) (4.5) 

On the other hand, using Liischer's results [37] it only differs from the infinite volume mass 
in terms exponentially falling with L. Using eqn. (4.3) 

Am(L) = -F^ (iTT, 0) + 0(6-"^) 
m 
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Similarly, the vacuum expectation value receives only corrections falling off exponentially with 
L. Therefore we obtain 

mmmL = ^ (F*(m,o) + mmm) +... 

with the ellipsis denoting residual finite size corrections. Note that the factor mL is just the 
one-particle Bethe-Yang Jacobian pi{9) = mLcosh^ evaluated for a stationary particle = 0. 

We can extend the above result to moving particles in the following way. Up to residual 
finite size corrections, the one-particle energy is given by 

E{L) = + 

with 

2tts 

where s is the Lorentz spin (which is identical to the particle momentum quantum number). 
Therefore 

E^E = mAm 

whereas perturbation theory gives: 

AE = XLi{{sm{s})L-{0\mL) 

and so we obtain 

(WI^IW)l = ^ (F''{i7^,0)+p^{emm) +■■■ (4.6) 
Pi{9) V / 

where ^ 

sinh 6 = — - =^ pi (9) = y/mPlX+47r^ 
mL 

Figure (4.1) shows the comparison of eqn. (4.6) to numerical data obtained from Lee- Yang 
TCSA: the matching is spectacular, especially in the so-called scaling region (the volume range 
where residual finite size corrections are of the order of truncation errors, cf. [16]) where the 
relative deviation is less than 10~^. Here and in all following plots we use the dimensionless 
volume parameter I = mL, and the matrix elements are also measured in units of m (cf. [16] 
for details). Diagonal one-particle matrix elements for the Ising model are shown in figure 4.2, 
where we similarly use natural units given by the mass m = m\ of the lightest particle Ai, 
just as in all subsequent plots related to the Ising model. 

One can use a similar argument to evaluate diagonal two-particle matrix elements in finite 
volume. Let us assume that the theory considered has diagonal scattering as in section 2.1. 
The two-particle Bethe-Yang equations remain valid even in a non-integrable theory as long 
as the total energy of the two-particle state remains under the inelastic threshold [38], and 
therefore the energy levels can be calculated from 

mil L sinh e'l + K^i - ^2) = 27r/i 
rui^ L sinh 6*2 + 5(6*2 - 6*1 ) = 27rJ2 

and (up to residual finite size corrections) 

E2{L) = E2pt{L) - Eo{L) = rui^ cosh^i + rrii^ cosh 6*2 
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/ll 




I ^ I 

5 10 15 20 25 30 

Figure 4.1: Diagonal 1-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 



where ii and 12 label the particle species. After a somewhat tedious, but elementary calculation 
the variation of this energy difference with respect to A can be determined, using (4.3) and 
(4.4): 

AE2iL) = A ^ ^(f^ (62 + ITT, 01 +177,91,92) + mil L cosh ^1^2* (m,0)i2i2 

Pni, [01,02)^ ^ ^ '''''''' 

+mi^ L cosh 92F^ {177,^)1-^^^ 

where all quantities (such as Bethe-Yang rapidities 9i, masses rr?,,- and the two-particle state 
density P2) are in terms of the A = theory. This result expresses the fact that there are two 
sources for the variation of two-particle energy levels: one is the mass shift of the individual 
particles, and the second is due to the variation in the interaction. On the other hand, in 
analogy with (4.5) we have 

/^E2{L) = XL{i,i,{{h,l2}mh,l2})ni,,L - (0|*|0)l) 
and so we obtain the following relation: 

i,i,{{Ii,l2}mh,l2})ni„L = ^^^(fI (92 + i7r,9i+i7r, 91,92) 



+mi-^ L cosh 91F2 (ivr, 0)^242 

+mi^Lcosh92Fi{iTr,0)i,i^ + (0|*|0)) + . . . (4.7) 
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Figure 4.2: Diagonal 1-particle matrix elements in the Ising model. The discrete points 
correspond to the TCSA data, while the continuous line corresponds to the prediction from 
exact form factors. 
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Figure 4.3: Diagonal 2-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 

where the ellipsis again indicate residual finite size effects. The above argument is a gener- 
alization of the derivation of the mini-Hamiltonian coefficient C in Appendix C of [15]. This 
formula is tested against numerical data in the Lee- Yang model in figure 4.3, and the agree- 
ment is as precise as it was for the one-particle case. Similar results can be found in the Ising 
case; they are shown in figure 4.4. 

4.2 Generalization to higher number of particles 

Let us now introduce some more convenient notations. Given a state 

\{h ■ ■■In})ii...in 

we denote 

p{{ki ,...,kr})L = P^,^ ...i,^ [K ,---Ar) (4.8) 

where 9i,l = 1, . . . ,n are the solutions of the n-particle Bethe-Yang equations (2.9) at volume 
L with quantum numbers I\,...,In and . . . ,kr}.,L) is the r-particle Bethe-Yang Jacobi 

determinant (2.10) involving only the r-element subset 1 < k\ < ■ ■ ■ < kr < n oi the n 
particles, evaluated with rapidities ^fcu • • • > ^fc,-- Let us further denote 

J^{{ki,...,kr})L = Firi^ki, ■ ■ ■ ,dkr)iki-ikr 

where 

FiniOi, On)h...i„ = lim (0n + in + e, . . . , Oi + iiT + €, Oi, . . . , en)n...i„i„...h (4-9) 
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Figure 4.4: Diagonal 2-particle matrix elements in the Ising model. The discrete points 
correspond to the TCSA data, while the continuous line corresponds to the prediction from 
exact form factors. 

is the so-called symmetric evaluation of diagonal n-particle matrix elements, which we analyze 
more closely in the next subsection. Note that the exclusion property mentioned at the end 
of subsection 2.1 carries over to the symmetric evaluation too: (4.9) vanishes whenever the 
rapidities of two particles of the same species coincide. 

Based on the above results, we conjecture that the general rule for a diagonal matrix 
element takes the form of a sum over all bipartite divisions of the set of the n particles 
involved (including the trivial ones when A is the empty set or the complete set {1, . . . ,n}): 



ii...in,L 



p{{l, . . .,n})L 

HA)lp{{1, ...,n}\A)L + ©(e-'^^ 



(4.10) 



Ac{l,2,...n} 



This rule can be tested against matrix elements with n = 3 and n = 4 in the Lee- Yang model, 
which are displayed in figures 4.5 and 4.6, respectively. The agreement is excellent as before, 
with the relative deviation in the scaling region being of the order of 10~^. 

5 Diagonal matrix elements in terms of connected form factors 

In this section we discuss diagonal matrix elements in terms of connected form factors, and 
prove that a conjecture made by Saleur in [21] exactly coincides with our eqn. (4.10). To 
simplify notations we omit the particle species labels; they can be restored easily if needed. 

5.1 Relation between connected and symmetric matrix elements 

The purpose of this discussion is to give a treatment of the ambiguity inherent in diagonal 
matrix elements. Due to the existence of kinematical poles (2.6) the expression 
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Figure 4.5: Diagonal 3-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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Figure 4.6: Diagonal 4-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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which is relevant for diagonal multi-particle matrix elements, is not well-defined. Let us 
consider the regularized version 

F2n{01 + ITT + €1,62 + ZTT + 62, ...,6'„ + ZTT + 02-1^1) 

It was first observed in [36] that the singular parts of this expression drop when taking the 
limits Cj — > simultaneously; however, the end result depends on the direction of the limit, 
i.e. on the ratio of the e, parameters. The terms that are relevant in the limit can be written 
in the following general form: 

F2niOi + iTT + ei, 6*2 + ^7^ + 62, 6'n + ztt + en, On, -, 6*2, 6*1) = (5.1) 

n ^ n n n 

JJ — • X X ••■ X "n«2...in(^i) • • • ,^n)enei2...ej„ + ... 

i=l * il=li2=l jn,=l 

where aiii2...i„ is a completely symmetric tensor of rank n and the ellipsis denote terms that 
vanish when taking ej — ^ simultaneously. 

In our previous considerations we used the symmetric limit, which is defined by taking all 
Ci equal: 

^2, -, On) = lim F2ni0i + iTr + e,02 + iTT + e, ^„ + + e, 0n, -, 02, 0i) 

It is symmetric in all the variables O\,...,0n- There is another evaluation with this symmetry 
property, namely the so-called connected form factor, which is defined as the independent 
part of eqn. (5.1), i.e. the part which does not diverge whenever any of the ej is taken to zero: 

F^^{0i,e2,...,en) = n\ai2...n (5.2) 
where the appearance of the factor n! is simply due to the permutations of the e^. 

5.1.1 The relation for n < 3 

We now spell out the relation between the symmetric and connected evaluations for n = 1, 2 
and 3. 

The n = 1 case is simple, since the two-particle form factor ^2(^1, ^2) has no singularities 
at Oi = 62 + ITT and therefore 

Fi{0) = = F2(iTr,0) (5.3) 

It is independent of the rapidities and will be denoted in the sequel. 
For n = 2 we need to consider 

7-. //I , • , /I , ■ , n n\ ailCi + 2ai2eie2 + 02262 /r ^\ 

^4(^1 + ZTT + €1,02 + nr + 62, 02, 0i) ~ (5.4) 

€162 

which gives 

Fi{9i,02) = aii + 2ai2 + a22 
FiiOi,02) = 2ai2 
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The terms an and 022 can be expressed using the two-particle form factor. Taking an in- 
finitesimal, but fixed £2 7^ 

Res F4(6i +177 + €1,02 + i-!r + €2, 02, Oi) = 02262 

€1=0 

whereas according to (2.7) 

Res F^Oi+iTT+ei, 02+1^+62, 62, 6*1) = i (1 - - ^2)5(^1 - ^2 - ^tt - €2)) F2(02+i7r+e2, ^2) 

€1=0 

To first order in £2 

s{ei -92-m- €2) = 8(62 -01 + €2) = s{e2 - ei)ii + 1^(62 - ^1)62 + . . . ) 

where 

is the derivative of the two-particle phase shift defined before. Therefore we obtain 

022 = ^{02 - ei)F^ 

and similarly 

an = <p{ei - e2)Fi 

and so 

Fi{ei,e2) = Fi{6i,e2) + 2^{ei - e2)F2{m, 0) (5.5) 

In the case of the trace of the energy-momentum tensor O the following expressions are known 
[24] 

Ff = 2TTm^ 

= 87rmV(^i - ^2) cosh^ 

F^'" = 47rmV(^i - ^2) cosh(0i - ^2) 

and they are in agreement with (5.5). 

For n = 3, a procedure similar to the above gives the following relation: 

Fi{ei,e2, Os) = F^iOi, 62, e^) + [fhoi, e2){^{ei - e^) + ip{e2 - e^)) + permutations] 

+3F| [ip{ei - e2)^{0i - O3) + permutations] (5.6) 
where we omitted terms that only differ by permutation of the particles. 

5.1.2 Relation between the connected and symmetric evaluation in the general 
case 

Our goal is to compute the general expression 

F2niOl, ■ ■ ■,On\f-l, . . . ,e„) = F2n{0i + ivT + €1,62 + 177 + 62, ...,^„ + 177 + en,On, ...,02,Oi) (5.7) 

Let us take n vertices labeled by the numbers 1,2, ... ,n and let G be the set of the directed 
graphs Gi with the following properties: 

• Gi is tree-like. 

• For each vertex there is at most one outgoing edge. 
For an edge going from i to j we use the notation Eij. 
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Theorem 1 (5.7) can be evaluated as a sum over all graphs in G, where the contribution of 
a graph Gj is given by the following two rules: 

• Let Ai = {ai,a2, . . . ,a^} be the set of vertices from which there are no outgoing edges 
in Gj. The form factor associated to Gj is 

-^2m(^ai ! ^'a2 ' ■ ■ ■ ' ^am) (5-8) 

• For each edge Ej^ the form factor above has to be multiplied by 

Note that since cannot contain cycles, the product of the ei/ej factors will never be trivial 
(except for the empty graph with no edges). 

Proof The proof goes by induction in n. For n = 1 we have 

Fi{ei) = Fi{e^) = F2{m,^) 

This is in accordance with the theorem, because for n = 1 there is only the trivial graph which 
contains no edges and a single node. 

Now assume that the theorem is true for n — 1 and let us take the case of n particles. 
Consider the residue of the matrix element (5.7) at e„ = while keeping all the finite 

R = Res F2n{0i..en\f^l--^n) 

en=0 

According to the theorem the graphs contributing to this residue are exactly those for which 
the vertex n has an outgoing edge and no incoming edges. Let Rj be sum of the diagrams 
where the outgoing edge is Enj for some j = l,...,n — 1, and so 

n— 1 

R = Y,Rj 

J=l 

The form factors appearing in Rj do not depend on On- Therefore we get exactly the diagrams 
that are needed to evahiate -F2(n-i)(^i"^n-i|ei"en-i), apart from the proportionality factor 
associated to the link Enj and so 

Rj = —(p{Oj — ^„)F2(„_i)(6'i..^„_i|ei..e„_i) 

and summing over j gives 
R = {eMOl - On) + €2(^(^2 -On) + --- + e„_i^(0„_i - ^n))i^2(n-l) (^1 --^n-l kl-en-i) (5.9) 

In order to prove the theorem, we only need to show that the residue indeed takes this form. 
On the other hand, the kinematical residue axiom (2.6) gives 

= ^ f 1 - n - ^o)S{On - Oj -m- e,) j F2(n-l)(^l-^n-l|ei..en-l) 
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Figure 5.1: The graphs relevant for n = 2 
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Figure 5.2: The graphs relevant for n = 3 

which is exactly the same as eqn. (5.9) when expanded to first order in ej. 

We thus checked that the theorem gives the correct result for the t erms that include a 1/ e^j 
singularity. Using symmetry in the rapidity variables this is true for all the terms that include 
at least one 1/ej for an arbitrary i. There is only one diagram that cannot be generated by 
the inductive procedure, namely the empty graph. However, there are no singularities (l/e^ 
factors) associated to it, and it gives -F2n(^i' • • • ' ^n) by definition. Qed. 

We now illustrate how the theorem works. For n = 2, there are only three graphs, depicted 
in figure 5.1. Applying the rules yields 

F4{ei,92\ei,e2) = Fi{9i, 82) + (^(^1 - ^2) (- + - 

which gives back (5.5) upon putting ei = €2. For n = 3 there are 4 different kinds of graphs, 
the representatives of which are shown in figure 5.2; all other graphs can be obtained by 
permuting the node labels 1,2,3. The contributions of these graphs are 



(a) : F^{e,,e2,es) 

(6) : -¥^(^1- 02)^1(^2, ^3) 

(c) : --ip{9i - 92Me2 - es)F^ = %{ei - e2M92 - es)Fi 

ei 62 ei 

(d) : --^i9i - e2Me3 - e2)F^ 

Adding up all the contributions and putting ei = £2 = £3 we recover eqn. (5.6). 
5.2 Consistency with Saleur's proposal 

Saleur proposed an expression for diagonal matrix elements in terms of connected form factors 
in [21], which is partially based on earlier work by Balog [39] and also on the determinant 
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formula for normalization of states in the framework of algebraic Bethe Ansatz, derived by 
Gaudin, and also by Korepin (see [40] and references therein). To describe it, we must extend 
the normalization of finite volume states defined in [16] to the case when the particle rapidities 
form a proper subset of some multi-particle Bethe- Yang solution. 

According to [16], the normalization of a finite volume state is given by 

\{h, ■ ■ ■ ,In])L = I 1^1, ■ ■ ■ ,Gn) 

in terms of the infinite volume state with rapidities 6i,. . . ,6n, which are the solutions of the 
Bethe- Yang equations (2.9) for the given quantum numbers /i, . . . ,/„ at volume L (we again 
omit the particle species labels, and also denote the n-particle determinant by p„). Let us 
take a subset of particle indices A E {1, . . . ,n} and define the corresponding sub-determinant 
by 

Pn{ei,...,en\A) = dctjjl''^ 

where iT^"'' is the sub-matrix of the matrix J'^"'^ defined in eqn. (2.10) which is given by 
choosing the elements whose indices belong to A. The full matrix can be written explicitly as 



j{n) ^ 



( EiL + ipi2-\ \-(pin -^12 ■ ■ ■ -V^ln 

-¥'12 E2L 9?21 + ¥'23 H ^ ¥'2n • • • -f2n 

\ — V'ln —V2n ■■■ EnL + ipin-\ h^n-l^nj 



where the following abbreviations were used: Ei = mjcosh^j, ipij = ipji = ip{Oi — 6j). Note 
that pn depends on all the rapidities, not just those which correspond to elements of A. It is 
obvious that 

Pn{dl, ...,9n) = Pn{dl, ■ ■ ■ , On\{^, •••,"-}) 

Saleur proposed the definition 

{{Ok}keA\{ek}keA)L = PniOu ■ ■ ■ , en\A) (5.10) 

where 

\{(^k}keA)L 

is a "partial state" which contains only the particles with index in A, but with rapidities that 
solve the Bethe- Yang equations for the full n-particle state. Note that this is not a proper 
state in the sense that it is not an eigenstate of the Hamiltonian since the particle rapidities do 
not solve the Bethe- Yang equations relevant for a state consisting of particles (where |^| 
denotes the cardinal number - i.e. number of elements - of the set ^4). The idea behind this 
proposal is that the density of these partial states in rapidity space depends on the presence 
of the other particles which are not included, and indeed it is easy to see that it is given by 

Pn(^l,---,^nU). 

In terms of the above definitions, Saleur's conjecture for the diagonal matrix element is 

1 

E F^^^0k}keA)m, ■ ■ -JnlA) + 0(6-^^^) 

Ac{l,2,...n} 



n...iA{h...Inm{h...In})i,...ir.,L = (5-11) 
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which is just the standard representation of the full matrix element as the sum of all the 
connected contributions provided we accept eqn. (5.10). The full amplitude is obtained by 
summing over all possible bipartite divisions of the particles, where the division is into particles 
that are connected to the local operator, giving the connected form factor F'^ and into those 
that simply go directly from the initial to the final state which contribute the norm of the 
corresponding partial multi-particle state. 

Using the results of subsection 5.1, it is easy to check explicitly (which we did up to n = 3) 
that our rule for the diagonal matrix elements as given in eqn. (4.10) is equivalent to eqn. 
(5.11). We now give a complete proof for the general case. 

Theorem 2 

X] ^2\A\i{^k}keA)p{^l, ■ ■ ■ ,dn\A) = X F2\A\i{^k}keA)p{{dk}keN\A) (5-12) 
AcN AcN 

where we denoted TV = {1, 2, ... , n}. 

Proof The two sides of eqn. (5.12) differ in two ways: 

• The form factors on the right hand side are evaluated according to the „symmetric" 
prescription, and in addition to the connected part also they contain extra terms, which 
are proportional to connected form factors with fewer particles. 

• The densities p on the left hand side are not determinants of the form (2.10) written 
down in terms of the particles contained in \ ^: they contain additional terms due to 
the presence of the particles in A as well. 

Here we show that eqn. (5.12) is merely a reorganization of these terms. 

For simplicity consider first the term on the left hand side which corresponds to A = 
{m + 1, m + 2, . . . , 77,}, i.e. 

^2mi^m+l, ■ ■ ■ , dn)p{^l, • • • , ^n|^) 

We expand p in terms of the physical multi-particle densities p. In order to accomplish this, 
it is useful to rewrite the sub-matrix 4 as 



\ i=m+l J 

where is the m-particle Jacobian matrix which does not contain any terms depending on 
the particles in A. The determinant of J^\^j^ can be written as a sum over the subsets of 



i=m+l 



i=m+l 
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N\A. For a general subset B C iV \ ^ let us use the notation B = {61, 62, • • • , We can 

then write 

p(ei,...,e„|^) = detj(«)U\A = E 

B 

where p{N \ U B)) is the /O-density (2.10) written down with the particles in N \ {AL) B). 

Applying a suitable permutation of variables we can generalize eqn. (5.13) to an arbitrary 
subset Ac N: 

\B\ 

~p{ei, er,\A) = det \n\a = 5] p(iV \ (^ U B)) Y^iU <Pb„cJ (5.14) 

B C i=l 

where the second summation goes over all the sets C = {ci,C2, . . . )C|b|} with \C\ = \B\ and 
Ci G A. The left hand side of eqn. (5.12) can thus be written as 

^ p{N\{AUB))J2F{A,b,C){5-'^5) 
A,BcN ^ 

\B\ 

F{A,B,C) = P2\A\i{^k}k&A)Y[y:>bi,Ci 

1=1 

We now show that there is a one-to-one correspondence between all the terms in (5.15) and 
those on the right hand side of (5.12) if the symmetric evaluations are expanded according 
to Theorem 1. To each triplet {A,B,C) let us assign the graph G(a,b,c) defined as follows: 

• The vertices of the graph are the elements of the set AU B. 

• There are exactly \B\ edges in the graph, which start at 6j and end at Cj with i = 
1,...,\B\. 

The contribution of G(yi,_B,c) to F2Qj^^_^^^^-^{{6k}keAuB) is nothing else than i^(^,_B,c) which 
can be proved by applying the rules of Theorem 1. Note that all the possible diagrams with 
at most n vertices are contained in the above list of the G(yi,_B,c) 1 because a general graph G 
satisfying the conditions in Theorem 1 can be characterized by writing down the set of vertices 
with and without outgoing edges (in this case B and A) and the endpoints of the edges (in 
this case C). 

It is easy to see that the factors p{N \ (^ U B)) multiplying the iii (5.15) are also 

the correct ones: they are just the density factors multiplying F^^^^^^^^^-j{{Ok}keAuB) on the 
right hand side of (5.12). Qed. 

6 Zero-momentum particles 
6.1 Scaling Lee- Yang model 

In the scaling Lee- Yang model, with a single type of particle, there can only be a single particle 
of zero momentum in a multi-particle state due to the exclusion principle. For the momentum 



piN\{AuB))l[ 



n 

E 

i=m+l 



(5.13) 



E^2|A|(WfeeA)^(^l>--->^n|^) = 

AcN 



where 
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to be exactly zero in finite volume it is necessary that the all other particles should come with 
quantum numbers in pairs of opposite sign, which means that the state must have 2n + 1 
particles in a configuration 

\{h, ■■■,In,^, —In, ■■■ 1 —h})L 
Therefore we consider matrix elements of the form 

{{![,..., 4, 0, -4, . . . , -limih, ...,/;, 0, -ii, -h})L 

(with A; = or Z = corresponding to a state containing a single stationary particle). W'e also 
suppose that the two sets {/i, . . . ,1^} and . . . are not identical, otherwise we have 
the case of diagonal matrix elements treated in section 4. 
We need to examine form factors of the form 

i^2fc+2i+2(«7r + ...,i'K + e'^,m-0'j,,...,m- e[,m + e,Q,0u . . .,9i,-ei, . . . 

where the particular ordering of the rapidities was chosen to ensure that no additional S matrix 
factors appear in the disconnected terms of the crossing relation (2.2). Using the singularity 
axiom (2.6), plus unitarity and crossing symmetry of the S'-matrix it is easy to see that the 
residue of the above function at ^ = vanishes, and so it has a finite limit as 6^0. However, 
this limit depends on direction just as in the case of the diagonal matrix elements considered 
in section 4. Therefore we must specify the way it is taken, and just as previously we use 
a prescription that is maximally symmetric in all variables: we choose to shift all rapidities 
entering the left hand state with the same amount to define 

\imF2k+2i+2{iT^ + O'l + e, . . . ,iTT + e'l^ + e,iTr - + e, 

e— ♦O 

i7r + e,0, 91,..., 01,-01,..., -Oi) 

Using the above definition, by analogy to (4.10) we conjecture that 

f2k+i,2i+i = {{![,..., l'k,0,-li...,-l[m{h,...,Ii,0,-Ii,...,-h})L (6.2) 

1 ^ 

^P2k+i{e[, ~ei 0, -~e'^, -e[)p2i+i{9i, ~0u o, -~ei,..., -h) 
(j'kA^'i^ ■ ■ ■ ' ^ll^i> ■■■A) + mLF2k+2i{iTT + e[,...,m + e'k, 

i7r-9'^,...,i7r- e[,9i, . . . ,9i,-9i, . . . , -^i)) + 0(6-^^^) 

where 9 denote the solutions of the appropriate Bethe-Yang equations at volume L, /?„ is a 
shorthand notation for the n-particle Bethe-Yang density (2.10) and equality is understood up 
to phase factors. We recall from our previous work [16] that relative phases of multi-particle 
states are in general fixed differently in the form factor bootstrap and TCSA. Also note that 
reordering particles gives phase factors on the right hand side according to the exchange 
axiom (2.4). This issue is obviously absent in the case of diagonal matrix elements treated in 
sections 4 and 5, since any such phase factor cancels out between the state and its conjugate. 
Such phases do not affect correlation functions, or as a consequence, any physically relevant 
quantities since they can all be expressed in terms of correlators. 



. . . ,nr ■ 



9[ + e, 



(6.1) 



25 



There is some argument that can be given in support of eqn. (6.2). Note that the zero- 
momentum particle occurs in both the left and right states, which actually makes it unclear 
how to define a density similar to p in (5.10). Such a density would take into account the 
interaction with the other particles. However, the nonzero rapidities entering of the two states 
are different and therefore there is no straiglitforwaxd way to apply Saleur's recipe (5.11) here. 
Using the maximally symmetric definition (6.1) the shift e can be equally put on the right hand 
side rapidities as well, and therefore we expect that the density factor multiplying the term 
P2k+2i in (6.2) would be the one-particle state density in which none of the other rapidities 
appear, which is exactly mL for a stationary particle. This is a natural guess from eqn. 
(4.10) which states that when diagonal matrix elements are expressed using the symmetric 
evaluation, only densities of the type p appear. 

Another argument can be formulated using the observation that eqn. (6.2) is only valid if 
J^k,l is defined as in (6.1); all other possible ways to take the limit can be related in a simple 
way to this definition and so the rule (6.2) can be rewritten appropriately. Let us consider 
two other natural choices 

lim ^2fe+2i+2(«7r + e[, . . . ,iTr + e'k,iTr - e'f., . . . ,iTr - e[,iTr,e,ei, . . . ,61, -Ou -Oi) 

e— »0 

^i^iid'l, ■ ■ ■,dk\dl, ■■■,dl) = 

lini F2fc+2;-h2(«7r + e[, . . . ,177 + 9'k,iTr - 9'^, . . . ,117 - e[,m + e, 0, ^i, . . . , 9^ -9i, -9i) 

in which the shift is put only on the zero-momentum particle on the right/left, respectively. 
Using the kinematical residue axiom (2.6), can be related to via 

J^k,i{di, ■ ■ ■,0k\^i, = J^k,A9[, . . .,9'j^\9i, ...,9i) 

I 

+2Y,f{0i)F2k+2i{iT^ + . . . , ivr + ivr - . . . , ZTT - ^1, . . . , Ou -Ou • • • , -^i) 
j=i 

^k,l{^l, ■ ■ ■ j^fel^l) ■ ■ ■ = ^k^ii^l^ ■ ■ ■ J^jfcl^lj ■ ■ ■ i^l) 

k 

-2Y,^{9[)F^k+2i{iT^ + 9'^,...,m + 9lm-9l...,m- 9[, 9i, . . . A, -9u • • • , -^1) 

i=l 

With the help of the above relations eqn. (6.2) can also be rewritten in terms of The way 
T and therefore also eqn. (6.2) are expressed in terms of shows a remarkable and natural 
symmetry under the exchange of the left and right state (and correspondingly J^'^ with J^~), 
which provides a further support to our conjecture. 

The above two arguments cannot be considered as a proof; we do not have a proper 
derivation of relation (6.2) at the moment. On the other hand, as we now show it agrees very 
well with numerical data which would be impossible if there were some additional ip terms 
present; such terms, as shown in our previous work [16] would contribute corrections of order 
1/i in terms of the dimensionless volume parameter I = mL. 

Data for the case of 1-3 and 3-3 matrix elements are shown in figures 6.1 and 6.2, respec- 
tively. In order to strengthen the support for eqn. (6.2) we must find 5-particle states. This 
is not easy because they are high up in the spectrum, and identification using the process of 
matching against Bethe-Yang predictions (as described in [16]) becomes ambiguous. We could 
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Figure 6.1: l-particle-3-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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Figure 6.2: 3-particle-3-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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Figure 6.3: Identifying the 5-particle state using form factors. The discrete points correspond 
to the TCSA data, while the continuous Hue corresponds to the prediction from exact form 
factors. 



identify the first 5-particle state by combining the Bethe-Yang matching with predictions for 
matrix elements with no disconnected pieces given by eqn. (2.11), as shown in figure 6.3. Some 
care must be taken in choosing the other state because many choices give matrix elements that 
are too small to be measured reliably in TCSA: since vector components and TCSA matrices 
are mostly of order 1 or slightly less, getting a result of order 10"^ or smaller involves a lot of 
cancellation between a large number of individual contributions, which inevitably leads to the 
result being dominated by truncation errors. Despite these difficulties, combining Bethe-Yang 
level matching with form factor evaluation we could identify the first five-particle level up to 
I = 20. 

The simplest matrix element involving a five-particle state and zero-momentum discon- 
nected pieces is the 1-5 one, but the prediction of eqn. (6.2) turns out to be too small to be 
usefully compared to TCSA. However, it is possible to find 3-5 matrix elements that are suffi- 
ciently large, and the data shown in figure 6.4 confirm our conjecture with a relative precision 
of somewhat better than 10^"' in the scaling region. 

We close by noting that since the agreement is better than one part in 10'^ in the scaling 
region, which is typically found in the range of volume Z ~ 10 ... 20, and also this precision 
holds for quite a large number of independent matrix elements, the presence of additional 
terms in eqn. (6.2) can be confidently excluded. 

6.2 Ising model in magnetic field 

In figure 6.5 we show how the prediction (6.2) describes a 1-3 matrix element in the Ising 
model; since all particles in this example are of species Ai, the formula carries over without 
essential modifications. 

However, due to the fact that the Ising model has more than one particle species, it is 
possible to have more than one stationary particles in the same state. Our TCSA data allow 
us to locate one such state, with a stationary Ai and A2 particle, and extending our previous 
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Figure 6.4: 3-particle-5-particle matrix elements in the scaling Lee- Yang model. The discrete 
points correspond to the TCSA data, while the continuous line corresponds to the prediction 
from exact form factors. 
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Figure 6.5: Ai — AiAiAi matrix element in Ising model with a zero-momentum particle 
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Figure 6.6: Ai — A1A2 matrix element in Ising model with zero-momentum particle 



considerations we have the prediction 

/i,i2 = i({0}|*|{0,0})i2 = ^ \ ^ (limF3(i7r + e,0, 0)112 +miLFi (0)2) 

where ^1(0)2 is the one-particle form factor corresponding to A2- This is compared to TCSA 
data in figure 6.6 and a convincing agreement is found. 

Note that in both of figures 6.5 and 6.6 there is a point which obviously deviates from the 
prediction. This is a purely technical issue, and is due to the presence of a line crossing close 
to this particular value of the volume which makes the cutoff dependence more complicated 
and so slightly upsets the extrapolation in the cutoff. We also remark that we cannot check 
further matrix elements at the moment, because the appropriate form factor solutions have 
not yet been computed. 



7 Finite temperature correlators 

In this section we show how a systematical low-temperature expansion for correlation functions 
can be developed using the results presented so far. Finite temperature correlation functions 
have attracted quite a lot of interest recently. Leclair and Mussardo proposed an expansion 
for the one-point and two-point functions in terms of form factors dressed by appropriate 
occupation number factors containing the TEA pseudo-energy function [20], based on a quasi- 
particle description motivated by the thermodynamic Bethe Ansatz. As discussed in the 
introduction, their proposal for the two-point function was shown to be incorrect by Saleur 
[21]; on the other hand, he also gave a proof of the Leclair-Mussardo formula for one-point 
functions based on the conjecture formulated in eqn. (5.11), provided the operator considered 
is the density of some local conserved charge. Since we proved that our formula (4.10) for 
diagonal matrix elements is equivalent to Saleur's conjecture, our results in section 4 can be 
considered as a very convincing numerical evidence for the correctness of his argument. 

Another proposal for finite-temperature one-point functions was made by Delfino [23], who 
attempted to express them in terms of free-particle occupation numbers and the symmetric 
evaluation of diagonal matrix elements. It was shown by Mussardo that this proposal is not 
correct using a counter example where it disagreed with the Leclair-Mussardo expansion [24]. 
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Furthermore, Castro- Alvaredo and Fring also argued [25] that two-point functions cannot 
be obtained by a simple dressing procedure analogous to the Leclair-Mussardo expansion for 
one-point functions. They argued that one needs a more drastic change in the form factor 
program. 

All these issues are connected to the problem of finding a proper definition of the discon- 
nected pieces. Prom the crossing relation (2.2), these are infinite for the form factors defined 
in infinite volume, and subtraction of such infinities must be made with care in order to obtain 
the correct finite pieces. Because of the above difficulties there is also a development in the di- 
rection of finite temperature form factors (for a review cf. [41]); with further development, this 
other line of thought can also give a very useful formulation of finite temperature correlation 
functions. 

Here we use the idea that putting the system into a finite volume L provides a regularization 
for the form factors, which can even be considered physical since in the real world there are no 
infinite systems^. Our expressions for the finite volume form factors are valid up to exponential 
corrections in the volume, which makes it clear that performing the calculation in finite volume 
and then taking the limit L — > oo we should recover the proper finite temperature correlation 
function. Here we present the computation for the case of the one-point function up to the 
first three nontrivial orders; the calculation gets complicated for higher orders, but the recipe 
is straightforward. On general theoretical grounds, it is quite clear that our approach should 
also apply to the two-point function, or indeed to any multi-point correlator, but in order 
to keep the exposition short we do not go into these details here and leave them to future 
investigations. 



7.1 Leclair-Mussardo series expanded 

The finite temperature expectation value of a local operator O is defined by 



{O) 



R 



Tr (e" 



-RH 



O) 



Tr {e-RH) 



where R = 1/T is the temperature dependent extension of the Euclidean time direction used 
in thermal quantum field theory and H is the Hamiltonian. To keep the exposition simple we 
assume that the spectrum contains a single massive particle of mass m. Leclair and Mussardo 
proposed the following expression for the low temperature (T <^ m, or equivalently mR S> 1) 
expansion of the above one-point function: 



n=0 



1 1 

n! {2-Ky 



+ e- 



-^2n(^l' ■•■)^n) 



(7.1) 



where F|„ is the connected diagonal form factor defined in eqn. (5.2) and e{6) is the pseudo- 
energy function, which is the solution of the thermodynamic Bethe Ansatz equation 



/do' 
— ip{e - 6'')log(l + e" 
27r 



(7.2) 



^There is actually a little subtlety here, since we impose periodic boundary conditions which are also 
nonphysical, but we make use of the old intuition that nothing can actually depend very much on the choice 
of the boundary condition if the system is very large and has a finite correlation length (i.e. a mass gap). 
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The solution of this equation can be found by successive iteration, which results in 

2iT 2 7 27r 

/rift' rift" 
27r 27r 

Using this expression, it is easy to derive the following expansion from (7.1) 
//y~i\R /rnx i f rpc f —mRcosh.0 —2mRcosh.0\ 

= + i 2^^H' ) 

+ - / — — (F|(^i, 02) + 2$(0i - e2)F^) e-"^«^°«h0ig-miicosh02 

2 J 27r 27r 

+0 (7.4) 

where (O) denotes the zero-temperature vacuum expectation value. The above result can also 
be written in terms of the symmetric evaluation (4.9) as 

(C?)^ (C?) + J ^g— m-RcoshS g— 2m-Rcosh0^ _j_ 

1 r Mide2j^s^^^^^^^^_mR(cosh9,+coshe,) ^ q /g-3mi?^ ^75^ 

2 J 27r 27r 

where we used relations (5.3) and (5.5). 

For completeness we also quote Delfino's proposal: 

3— mil cosh 



\ \ ^ " g-mit COSH y. 



F|„(ei,...,0„) (7.6) 

which gives the following result when expanded to second order: 

//n\R /m\ I f TPS f^—mRcoshO „—2mRcosh6\ . 

{U}d = \0) + j —t'2 \^ - e j + 

1 r d0ide2ps^^^^^^^^.mRicoshe^+coshe2) + Q U-3mR\ (7 7) 

2 7 27r 27r 

Note that the two formulae coincide with each other to this order, which was already noted in 
[23]. However, this is not the case in the next order. Obtaining the third order correction from 
the Leclair-Mussardo expansion is a somewhat lengthy, but elementary computation, which 
results in 

1 f ^^^^s/fl Q a \ -mil(cosh6>i+cosh6»2+cosh6»3) 

6 7 27r 27r 27r ^' ^' ^' 

/^^^^^K^j^ ^^^g-mii(cosh 61+2 cosh 6*2) _|_ [ ^^^ps^-SmRcoshdi 
27r 27r J 27r 

_i J ^^FMOi - 02)e-"^«('^°«hei+2coshe2) (7 8) 
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where we used eqns. (5.3, 5.5, 5.6) to express the result in terms of the symmetric evaluation. 
On the other hand, expanding (7.6) results in 

1 f ^^^^^^psfQ^ ^^Ng-mR(cosh6>i+cosh6»2+cosh6»3) 

6 J 27r 27r 27r ^ ' ' 

fd6id92 sfa \^-mii(cosh0i+2cosh6»2) I / '^^l ps^-Smii cosh 6»i /r, (^^ 

It can be seen that the two proposals differ at this order (the last term of (7.8) is missing from 
(7.9)), which was already noted by Mussardo using a toy model in [24], but our computation 
here is model independent and shows the general form of the discrepancy. We also need the 
third order correction explicitly so that we can compare it to the result of the computation 
performed in the next section. 

7.2 Low-temperature expansion for one-point functions 

We now evaluate the finite temperature expectations value in a finite, but large volume L: 

where Hi is the finite volume Hamiltonian, and Til means that the trace is now taken over 
the finite volume Hilbert space. For later convenience we introduce a new notation: 

1^*1, • • • ,On)L = \{h, ■ ■ ■ ,In})L 

where 9i,...,9n solve the Bethe-Yang equations for n particles with quantum numbers /i, . . . , 
at the given volume L. We can develop the low temperature expansion of (7.10) in powers of 
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+0(e-^"*^) (7.11) 



em ^ 0(2)^^(2) 
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The denominator of (7.10) can then be easily expanded: 




The primes in the multi-particle sums serve as a reminder that there exist only states for 
which all quantum numbers are distinct. Since we assumed that there is a single particle 
species, this means that terms in which any two of the rapidities coincide are excluded. All 
n-particle terms in (7.11) and (7.12) have a 1/n! prefactor which takes into account that 
different ordering of the same rapidities give the same state; as the expansion contains only 
diagonal matrix elements, phases resulting from reordering the particles cancel. The upper 
indices of the rapidity variables indicate the number of particles in the original finite volume 
states; this is going to be handy when replacing the discrete sums with integrals since it keeps 
track of which multi-particle state density is relevant. 

We also need an extension of the finite volume matrix elements to rapidities that are 
not necessarily solutions of the appropriate Bethe-Yang equations. The required analytic 
continuation is simply given by eqn. (4.10) 



where we made explicit the volume dependence of the n-particle density factors. The last term 
serves as a reminder that this prescription only defines the form factor to all orders in 1/L 
(i.e. up to residual finite size corrections), but this is sufficient to perform the computations 
in the sequel. 

Using the leading behaviour of the n-particle state density, contributions from the n- 
particle sector scale as L", and for the series expansions (7.11), (7.12) and (7.13) it is necessary 
that mL <^ e"*^. However if mR is big enough there remains a large interval 



where the expansions are expected to be valid. After substituting these expansions into (7.10) 
we will find order by order that the leading term of the net result is 0{L^), and the corrections 
scale as negative powers of L. Therefore in (7.10) we can continue analytically to large L and 
take the L — ^ oo limit. 

7.2.1 Corrections of order e~"^^ 

Substituting the appropriate terms from (7.13) and (7.11) into (7.10) gives the result 




(7.14) 



1 < mL < e' 



(O)^ = {0)l + E 
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Taking the -L — ^ oo limit one can replace the summation with an integral over the states in 
the rapidity space: 

i 

and using (4.6) we can write 

mO\e)L - {0)l) = Fi + 0(6-^^^) (7.15) 

so we obtain 

J 27r 

which coincides with eqn. (7.5) to this order. 
7.2.2 Corrections of order e"^"*^ 

Substituting again the appropriate terms from (7.13) and (7.11) into (7.10) gives the result 
0(1) 

\ ( \ 

I VP I 



-3mR\ 



9(2) n(2) 



The 0(e ^mi?^ terms can be rearranged as follows. We add and subtract a term to remove 
the constraint from the two-particle sum: 



5)(2)^g(2) 



(2) (2) 

The 9\ = 62 terms correspond to insertion of some spurious two-particle states with equal 
Bethe quantum numbers for the two particles (/i = 12). The two-particle Bethe-Yang equa- 
tions in this case degenerates to the one-particle case (as discussed before, the matrix elements 
can be defined for these "states" without any problems since we have the analytic formula 
(7.14) valid to any order in 1/L). This also means that the density relevant to the diago- 
nal two-particle sum is pi and so for large L we can substitute the sums with the following 
integrals 

^(1) fl(2)-fl(2) ff(2) fl(2) 

"1^2 "1 ""2 "1 '"2 
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Let us express the finite volume matrix elements in terms of form factors using (4.6) and (4.7): 

P2(^l,02) ((ef\^?f )l - {0)l) 

-Pi {0i) Pi (02) {{0i\O\9i)l + {e2\0\92)L - 2{0)l) = Fl{9i,e2) + 0(6"'^^) 
Combining the above relation with (7.15), we also have 

{6, e\o\e, e)L - = + o(e-^^) 

where we used that Fl(6, 9) = 0, which is just the exclusion property mention after eqn. (4.9). 
Note that 

=1 + 0(L-) 



and therefore in the limit L — ^ oo we obtain 

_ I d0 ^-2mRcosh.e ps I ^ f '^^Ips^Q^ ^^■jg-m_R(cosh6li+coshe2) 

7 27r 2 7 27r 27r 

which is equal to the relevant contributions in the Leclair-Mussardo expansion (7.5). 

7.2.3 Corrections of order e~^"^^ 

This calculation is rather long, and so it is relegated to the appendix. The net result is 

1 /"^^^^s//) Q ^ \g-mil(cosh6>i+cosh6l2+cosh6»3) 

6] 27r 27r 27r ^' ^' ^' 

^ r '^^^^ps^Q^ ^^■jg-mi?(coshei+2cosh6»2) _|_ f ps^-3mRcoshei 

7 27r 27r J 2n 

_1 j ^^Fi<f{ei - ^2)e-'"«(-^^^i+2coshe2) ^7 

which agrees exactly with eqn. (7.8). 
7.3 Remarks 

There are a few remarks which we wish to make. First, we see that the proposals by Leclair 
and Mussardo and by Delfino differ at the order e~^™^. The reason for this difference can be 
understood in the formalism developed here. Namely, the expansions (7.11) and (7.13) both 
contain positive powers of L. On physical grounds, they are expected to cancel completely 
order by order in the expansion. However, the state densities p depend on the interaction 

as well. This dependence is of order L^"^ , and it actually characterizes the ambiguity in the 
definition of the diagonal matrix element resulting from the resolution of the singularity (see 
eqn. (5.1)). Naively it drops out in the L — > oo limit, but actually some of these terms is 
multiplied by a positive L power from (7.13). In our derivation we evaluated every relevant 
contribution to all orders in 1/L (i.e. we only neglected residual finite size corrections). As a 
result, we could take the limit L ^ oo properly and get the correct finite part of the resulting 
expression. 

Taking this line of thought further, note that the leading term of every multi-particle 
density (whether it is degenerate in the sense defined in the appendix, or not) is always a 
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product of EiL factors where i runs over the number of particles and Ei is their energy. 
Therefore density terms whose leading behaviour is do not contribute explicit if factors. 
As far as there are only contributions of this type, the expansion of the one-point function, 
when written in terms of is just the same as in a free field theory. Indeed in the free 
field limit the Leclair-Mussardo expansion and the Delfino proposal are identical, since the 
pseudo-energy function is just e(^) = mR cosh and = (more generally, due to the 
absence of kinematical singularities the — > limit of (5.1) is independent of the direction). 

To have terms that depend explicitly on the interaction we need density contributions that 
naively scale as a positive power of L. When combining all such terms at a given order, the 
leading term must drop out, and the final result can only have a behaviour at large L. It 
is clear from our calculation detail above and in the appendix that the first order at which 
such an anomalous contribution arises is that of e~^"*^. Up to that order every individual 
term is finite as L — ^ oo. However, at third order there appear some "anomalous" density 
terms, namely those collected in (A. 7), which individually grow linearly in L. As required by 
general principles, the linear contribution cancels between them and so the L ^ oo limit is 
well-defined. However, the subleading terms always contain dependence on (p, and indeed they 
all vanish for a free theory (when <^ = 0), therefore it is only such terms that can contribute 
explicit (p dependence in the expansion. As a result, there remains an "anomalous" term which 
is just (—1 times) the derivative of the phase shift, and leads to the correction (A. 8), which is 
exactly the term absent in Delfino's expression. 

Strictly speaking, the above discussion is only valid if the expansion is written in terms of 
the symmetric evaluation F|„ ; rewriting it in terms of the connected form factors obviously 
introduces further (p dependence. As shown in the above argument, the real difference between 
the free and the interacting case can be properly observed when the expansion is written in 
terms of F^^, therefore it seems a more natural choice than using the connected form factors, 
as the behaviour specific to interacting theories can be seen much more clearly. 

Another important point is that our results give an independent support for the Leclair- 
Mussardo expansion. It is known that it coincides precisely with the exact TEA result for the 
trace of the energy-momentum tensor [20], and Saleur presented an argument for its validity 
when the operator considered is the density of a local conserved charge [21]. These arguments 
work to all orders, but only for a restricted set of local operators. On the other hand, our 
calculation above is model independent, and although we only worked it out to order e"^'"^, 
we expect that it coincides with the Leclair-Mussardo expansion to all orders. For a complete 
proof we need a better understanding of its structure, which is out of the scope of the present 
work. 

Furthermore, our method has a straightforward extension to higher point correlation func- 
tions. For example, a two-point correlation function 

{Oi{x)02mL [e-RHL) 

can be expanded inserting two complete sets of states 

TvL (e-^^^Oi(x)O2(0)) = 5^e-^^"(^)(n|O(x)|m)L(m|O(0)|n)L (7.17) 

m,n 

Since we now have a complete description of finite volume matrix elements to all orders in 1/L, 
the above expression can be evaluated along the lines presented in subsection 7.2, provided 
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that the intermediate state sums are properly truncated. We leave the explicit evaluation of 
expansion (7.17) to further investigations. 

Finally note that besides giving a systematic expansion in powers of e""'^, our method 
also gives the L dependence to all orders in 1/L (i.e. up to residual finite size effects), therefore 
it can also be used to study finite size corrections of correlators in the low temperature regime. 

8 Conclusions 

In this work we completed the description of finite volume matrix elements of local operators 
by considering those with disconnected pieces. There are two types of such matrix elements, 
namely (1) diagonal ones and (2) ones involving parity-invariant zero-spin states with zero- 
momentum particles. Our description is valid to any order in 1/L i.e. up to residual finite size 
corrections decaying exponentially with the volume L. The precise statements were formulated 
in subsection 2.3 and we then gave extensive numerical evidence for them. We also formulated 
and proved a general theorem relating the different possible evaluations of diagonal matrix 
elements, and showed that our results coincide with the proposal made by Saleur [21]. 

We then showed how to perform an expansion for finite temperature correlation functions, 
using the fact that finite volume acts as a regulator for the otherwise infinite disconnected 
pieces. The case we considered explicitly was that of one-point functions at finite temperature. 
We evaluated the first few orders in the low temperature expansion and showed that they 
coincide with the result conjectured by Leclair and Mussardo [20], but are different from 
Delfino's proposal [23] at third order. Some important aspects of this expansion were already 
discussed in subsection 7.3, which we do not repeat here. 

There is a number of interesting issues remaining. Our approach gives the finite volume 
form factors up to residual finite size effects, but combined with truncated conformal space 
one can achieve a precision of order 10~^ in the scaling Lee- Yang model, and 10~^ in the Ising 
model with magnetic field. It would be interesting to see how these results can be related 
to other approaches to finite volume form factors (see [42]) and whether the picture can be 
completed to give some sort of exact description in the case of integrable field theories. It also 
seems worthwhile to formulate a higher dimensional generalization of these results extending 
the approach of Lellouch and Liischer [43], which is expected to be relevant for lattice field 
theory. 

Another open issue is to give a more concise formulation of the finite temperature expansion 
discussed in section 7 that would make possible a partial resummation to recover the Leclair- 
Mussardo expression (7.1) which involves dressed form factors. 

It is even more interesting to write down the expansion for two-point correlators following 
the ideas outlined in subsection 7.3; a better method of organizing the contributions could 
be of great help here as well. Results for the two-point function can be compared e.g. to 
evaluation of correlation functions from truncated conformal space, and can also be used in 
further development of the finite temperature form factor program [41]. 
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A e ^"^^ corrections to the finite temperature one-point func- 
tion 

In order to shorten the presentation, we introduce some further convenient notations: 

Ei = m cosh 6i 

{eu...,en\o\ei,...,6n)L = {i...n\o\i...n)L 

Pn{9l,...,0n) = p{l...n) 

Summations will be shortened to 

E - E 

01.. Mn 1-n 

E' - E' 

6>i...6»„ l...n 

Given these notations, we now multiply (7.11) with (7.13) and collect the third order correction 
terms: 



+ 



1 ^ 'e-«(^i+-^2+^3) ((123|0|123)l - {0)l) 

^ 123 

{r^e-^^^^ ^Y^'^-n(E,+E,) ((23|0|23)i - (O)^) 



To keep trace of the state densities, we avoid combining rapidity sums. Now we replace the 
constrained summations by free sums with the diagonal contributions subtracted: 



E' 

12 


= E-: 

12 ] 


E 

L=2 


E' 

123 


= E- 

123 


\ 1=2,3 



1=2=3 

where the diagonal contributions are labeled to show which diagonal it sums over, but other- 
wise the given sum is free, e.g. 

E 

1=2,3 

shows a summation over all triplets , , where = 9^^ and ^3^^ runs free (it can 
also be equal with the other two). We also make use of the notation 

F{12...n) = Fi^i9i,...,9n) 
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so the necessary matrix elements can be written in the form 

p(123) ((123101123)2: - = F(123) + /9(1)F(23) + • • • + /9(12)F(3) + . . . 

p(122) ((122|C»|122)i - {0)l) = 2p{2)F{12) + 2p(12)F(3) + p(22)F(l) 

p(lll)((llip|lll)^-(0)i) = 3p(lll)F(l) 

p{12){{12\0\\2)l-{0)l) = F{12)+p{l)F{2)+p{2)F{l) 

p(11)((11|0|11U-(0)l) = 2p{l)F{l) 

p{1){{1\0\1)l-{0)l) = F{1) (A.l) 

where we used that F and p are entirely symmetric in all their arguments, and the ellipsis in 
the the first Une denote two phis two terms of the same form, but with different partitioning 
of the rapidities, which can be obtained by cyclic permutation from those displayed. We also 
used the exclusion property mentioned after eqn. (4.9). 

We can now proceed by collecting terms according to the number of free rapidity variables. 
The terms containing threefold summation are 



1 ^g-ii(Ei+B.+E3) ((i23p|123)i - {0)l) - ^5]5]((23p|23)i - {0)l) 

123 1 2,3 



^ 1,2 3 



Replacing the sums with integrals 



and using (A.l) we get 

b J ZTT ZTT ZTT 

+ / ^ff ^^"""^'^""^^'^^ (P(l)p(2)i^(3) - ^P(12)F(3)) 

where we reshuffled some of the integration variables. Note that all terms cancel except the 
one containing F(123) and writing it back to its usual form we obtain 

1 r ^^^^^^psf^Q^ 0^ ^^>jg-mR(cosh6>i+cosh6»2+cosh6>3) ^^2) 

6 7 27r 27r 27r ' ' 
It is also easy to deal with terms containing a single integral. The only term of this form is 

1 y2 e-^(-^i+-^2+-^3)({123|0|123)L-(0)L) 

1=2=3 
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(3) (3) (3) 

When all rapidities 61 ,02 ^ ^3 equal, the three-particle Bethe-Yang equations reduce to 
the one-particle case 

Therefore the relevant state density is that of the one-particle state: 



(3) 

mX sinh 9[ = 27r/i 



f ddl SrnR cosh ei ps /a n\ 

^ y (A.3j 

where we used that 
when L — > 00. 

The calculation of double integral terms is much more involved. We need to consider 
f E + E + E I e-«(^i+^^+^3) ((i23|0|123)i - {0)l) 

yi=2,3 1=3,2 2=3,1 J 

+^ E E e-^(^i+-^2+^3) ((23|C|23) l - {O) l) 

1 2=3 

+^ E E 6-^(^1+^^+^^) ((3|0|3)l - {0)l) (A.4) 
1=2 3 

We need the density of partially degenerate three-particle states. The relevant Bethe-Yang 
equations are 

ml sinh Oi + 5{ei -62) = 2itIi 
mL sinh 6*2 25(6*2 - ) = ^nh 

where we supposed that the first and the third particles are degenerate (i.e. Is = /i), and 
used a convention for the phase-shift and the quantum numbers where 6{0) = 0. The density 
of these degenerate states is then given by 



p{lS, 2) = det 



LEi+ ^{61-62) -^{01-02) 
-2^{0i - 02) LE2 + 2^{0i - 02) 



where we used that ip{0) = (p{~9). Using the above result and substituting integrals for the 
sums, we can rewrite eqn. (A.4) in the form 



+2 



+ 



1 1 '§^'^<<-'^''''''''''^^iMi)F(12)+2p(imi)+p{n)F(2)) 
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where the elhpsis denote two terms that can be obtained by cychcal permutation of the indices 
1,2,3 from the one that is exphcitly displayed, and these three contributions can be shown to 
be equal to each other by relabeling the integration variables: 

-\ I '^l^''"^'"'"^'^'''^^ (20(1)F(12) + 2p(12)F(l) +p(ll)F(2)) 

We first evaluate the terms containing F(23) which results in 
using that 

We can now treat the terms containing the amplitude F{\) = F{2) = F{3) = F.2. Exchanging 
the variables 6i ^ 62 in the second line and redefining ^3 02 in the third line of eqn. (A. 5) 
results in 

The combination of the various densities in this expression requires special care. Prom the 
large L asymptotics 

p{i)^EiL , p{ij)^EiEjL^ , p{ijk) ^ EiEjEkL^ , p{13,2) ^ EiE2L^ 

it naively scales with L. However, it can be easily verified that the coefficient of the leading 
term, which is linear in L, is exactly zero. Without this, the large L limit would not make 
sense, so this is rather reassuring. We can then calculate the subleading term, which requires 
tedious but elementary manipulations. The end result turns out to be extremely simple 

so the contribution in the L — >■ 00 limit turns out to be just 

2 J 27r 27r 

Summing up the contributions (A. 2), (A. 3), (A. 6) and (A.8) we indeed obtain (7.16). 
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